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Abstract 

Data system development activities for the Sea-viewing Wide Field-of-view Sensor (SeaWiFS) must begin well 
before the scheduled 1994 launch. To assist in these activities, it is essential to develop a simulated SeaWiFS 
data set as soon as possible. Realism is of paramount importance in this data set, including SeaWiFS spectral 
bands, orbital and scanning characteristics, and known data structures. Development of the simulated data 
set can assist in identification of problem areas that can be addressed and solved before the actual data are 
received. This paper describes the creation of the first version of the simulated SeaWiFS data set. The data set 
includes the spectral band, orbital, and scanning characteristics of the SeaWiFS sensor and SeaStar spacecraft. 
The information is output in the data structure as it is stored onboard. Thus, it is a level-0 data set which 
can be taken from start to finish through a prototype data system. The data set is complete and correct at the 
time of printing, although the values in the telemetry fields are left blank. The structure of the telemetry fields, 
however, is incorporated. Also, no account for clouds has been included. However, this version facilitates early 
prototyping activities by the SeaWiFS data system, providing a realistic data set to assess performance. 


1. INTRODUCTION 

The Sea-viewing Wide Field-of-view Sensor (SeaWiFS) 
is an ocean color sensor scheduled for launch in 1994. A 
follow-on to the highly successful Coastal Zone Color Scan- 
ner (CZCS) mission (1978-1986), SeaWiFS is intended to 
provide improved estimates of chlorophyll concentrations 
in the world's oceans. 

The SeaWiFS mission will provide an unprecedented 
set of global ocean color data on a regular basis (every two 
days in the absence of clouds, about a week otherwise). 
The CZCS, a proof-of-concept mission, was able to produce 
a single global image in its nearly eight year lifetime. The 
volumes of data produced by an operational global ocean 
color sensor pose challenges for a data processing system 
that were not encountered by the CZCS system. Develop- 
ment and rigorous testing activities must begin well before 
launch to avoid the sorts of delays in data processing and 
distribution that plagued the CZCS mission. 

To facilitate the development of the data system for 
processing SeaWiFS data, it is essential to develop a sim- 
ulated SeaWiFS data set as soon as possible. This data 
set should incorporate the characteristics of the SeaW r iFS 
sensor, including spectral bands and orbital and scanning 
characteristics, and should do so in a realistic manner to 
allow proper preparation for actual data processing. While 
many of the data characteristics, contents, and anomalies 
may not be known until just before launch, development of 
the simulated data set can assist in identification of prob- 
lem areas that can be addressed and solved before a fuller 
understanding of the data is achieved. 

This paper describes the details of the creation of the 
first version of the simulated SeaWiFS data set. The data 
set includes the spectral band, orbital, and scanning char- 
acteristics of the SeaWiFS sensor and SeaStar spacecraft. 
The information is output in the data structure as it is 
stored onboard. Thus, it is a level-0 data set which can be 


taken from start-to-finish through a prototype data sys- 
tem. The data set is complete and correct at the time of 
printing, although, the values in the telemetry fields are 
left blank (the structure of the telemetry fields is incor- 
porated, however). Also, no account for clouds has been 
included. 

CZCS pigment data and a variety of radiative transfer 
models were used to create realistic SeaWiFS total radi- 
ances. These total radiances were Earth-located using a 
SeaWiFS orbital model to provide SeaW T iFS-like viewing 
and solar geometries. 

2. BACKGROUND 

The improvements of SeaWiFS over CZCS are largely 
expected due to the existence of four new spectral bands: 
one at 412 nm (to discern dissolved organic matter, or 
Gelbstoff), one at 490 nm (to enable better estimates of 
chlorophyll at high concentrations), and two in the near- 
infrared at 765 and 865 nm (to improve the determination 
of atmospheric aerosols). Thus, SeaWiFS contains eight 
spectral bands compared to the four used in CZCS pro- 
cessing (CZCS actually contained six bands, but tw r o w^ere 
not used). Table 1 shows a comparison of the SeaWiFS 
spectral bands to the CZCS bands. 

In addition to different spectral bands, SeaWiFS is also 
in a different orbit than the CZCS and has different viewing 
geometries. A comparison between SeaWiFS orbital and 
sensor characteristics and those of the CZCS is summarized 
in Table 2. 

Due to limited onboard data recorder space, SeaWiFS 
will be able to obtain global coverage only at reduced res- 
olution (every fourth pixel along-track and along-scan). 
These data are called Global Area Coverage (GAC) data. 
High resolution data may be recorded in the remaining 
data recorder area. There is estimated to be about 20 
minutes per day of high resolution Local Area Coverage 
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(LAC) data. The effort here is to create both types of 
data, in order to simulate a representative SeaWiFS data 
set. 

Table 1. SeaWiFS spectral bands and center 
wavelengths, are shown with those of CZCS, for 
comparison. Wavelengths (A) are in nm. 


SeaWiFS 

CZCS \ 

Band No. 

A 

Band No. 

A 

i 

412 

1 

443 I 

2 

443 

2 

520 

3 

490 

3 

550 

4 

510 

4 

670 

5 

555 

5 

750 1 

6 

670 

6 

13,500 2 

7 

765 



8 

865 




1 . Intended for surface vegetation analyses only (Williams 
et al. 1985). 

2. Data suspect after 1979 (Williams et al 1985). 


3. METHODS 

In order to maximize the usefulness and realism of the 
simulated SeaWiFS data set, CZCS pigment data were 
chosen as the data source. The method for simulating 
SeaWiFS data involves creating both LAC and GAC data. 
Details of the method are described in the following, but 
a general overview here is helpful. 

The basic equation of radiative transfer for ocean color 
assumes separability of the radiance contributions from at- 
mospheric and oceanic components: 

Lt( A) = L r (A) + T(X)L g (X) -F L a (X) + T(A)Lvv*(A) (1) 


solar and viewing geometries, a realistic atmosphere was 
added by including Rayleigh scattering, aerosol scattering 
and absorption, water vapor absorption, oxygen absorp- 
tion, and sun glint reflectance. Radiance contributions 
of each were summed, and diffusely transmitted water- 
leaving radiance was added to produce simulated total ra- 
diance received by the sensor. Finally, data structures for 
telemetry were inserted and science data were formatted 
into 10-bit words to create level-0 data. 


3.1 Normalized Water-Leaving Radiances 

Creating normalized water-leaving radiances Lwn is 
the first step, because Lwn is independent of viewing and 
solar geometry. The method for creating Lwn is taken 
from Gordon et al. (1988a). According to this model, 


Lwn( A) 


(l-Pn)(l-PN)Fo(\)R{Q-,\) 

n 2 Q[l-rR(0~,\)} 


( 2 ) 


where p n is the Fresnel reflectance of the sea surface for 
normal incidence, p^ is a normalized mean value of surface 
reflectance for direct and diffuse irradiance for a flat sea, 
F 0 is the extraterrestrial irradiance corrected for Earth- 
sun distance, fl(0“,A) is the irradiance reflectance just 
below the sea surface, n is the index of refraction, Q is 
the irradiance-to- radiance ratio (equals tt for totally dif- 
fuse radiance), and r is the water- air reflectance for totally 
diffuse irradiance. 

From Fresnel’s Law, p n is known a priori (Jerlov 1976), 
Pn is estimated from Gordon et al. (1988a) at 0.043, n is 
taken to be 1.341 (Austin 1974), and r is taken to be 0.48 
(Gordon et al. 1988a). The ratio R{X)/Q is wavelength 
dependent, and is taken from Gordon et al. (1988a) as: 


where L t is the total radiance at the sensor, L r is the 
Rayleigh radiance, L g is the sun glint radiance, L a is the 
aerosol radiance, Lw is the water-leaving radiance, and 
T is the total transmittance (direct plus diffuse) from the 
ocean through the atmosphere to the spacecraft. The vari- 
able of interest in the creation of simulated data is L t , 
which is wdiat the sensor actually detects. To obtain re- 
alistic values of L t , with ranges and means observed in 


RW 

Q 



( 3 ) 


where &*>( A) is the spectral backscattering coefficient, and 
if (A) is the spectral attenuation coefficient, if (A) is deter- 
mined from Baker and Smith (1982) as 


if (A) = K w (\)+K c (\) + K g (\) (4) 


the real ocean, it is desirable to create it using represen- 
tative values of the components. This process begins with 
the generation of Lw derived from CZCS pigment values, 
and adding atmospheric contributions from simulated Sea- 
WiFS viewing and solar geometries. 

For both LAC and GAC data, real CZCS images of 
pigments were selected. The pigment fields were used to 
create normalized water-leaving radiances, Lwn{ A), at the 
eight SeaWiFS bands using a semi-analytic radiance model 
modified for SeaWiFS (Gordon et al. 1988a). Then, simu- 
lated SeaWiFS orbits were propagated through these sim- 
ulated radiance fields to produce actual SeaWiFS view- 


where ifvv(A) is the attenuation coefficient of pure seawa- 
ter, if c (A) is that for phytoplankton, and K g {\) is that for 
Gelbstoff, or yellow substances. 

Use of this model for a SeaWiFS simulation requires 
values for wavelength-dependent parameters be chosen for 
SeaWiFS bands. Values for the mean extraterrestrial ir- 
radiance F 0 were taken from Neckel and Labs (1984), and 
spectrally weighted over all of the SeaWiFS bandwidths 
assuming a full-width half-maximum (FWHM), Gaussian 
spectral response (Table 3). The values were corrected for 
Earth-sun distance following Gordon et al. (1983): 


ing geometries. An arbitrary day was chosen (the vernal _ f 2tt{D — 3)1 2 

equinox) to produce solar geometry. Given these simulated ^o(A) — F 0 ( A) 14- 0.0167 cos — - (5) 
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Table 2. Orbit and sensor characteristics for the SeaWiFS and the CZCS, shown for comparison. 


Characteristic 

CZCS 

SeaWiFS 

Altitude (km) 

955 

705 

Period (minutes) 

104.0 

98.9 

Inclination 

99.28° 

98.25° 

Equator Crossing Time (local) 

Noon 

Noon 

Node Type 

Ascending 

Descending 

Scan Width 

39.34° 

58.3° (LAG): 45° (GAG) 

Instantaneous Field of View (IFOV) 

0.05° 

0.09° 

Ground IFOV at Nadir (km) 

0.825 

1.12 

Pixels Along Scan 

1,968 

1,285 (LAC); 248 (GAC) 

Scan Period (seconds) 

0.134 

0.667 (LAC); 1.5 (GAC) 

Scan Plane Tilt 

±20° in 2° steps 

+20°. 0°, -20° 

Scan Ground Coverage (km) 

1,566 

2,802 (LAC); 1.502 (GAC) 

Maximum Spacecraft Zenith Angle 

46.8° 

70.8° (LAC); 51.7° (GAC) 

Digitization (bits) 

8 

10 


where Fo(A) is the mean extraterrestrial irradiance and D 
is the sequential day of the year. 

Table 3. Values for mean extraterrestrial irradi- 


ance Fo(A) (mW cm -2 ^m -1 ), k c ( A), ^(A), bw(X) 
(m -1 ), A(A) (m -1 ), and B( A) used in the creation 
of the SeaWiFS simulated data set. 


A 

F 0 

kc 

k' c 

bw 

A 

B 

410 

171.92 

0.208 

1.077 

0.0067 

0.00313 

0.20401 

443 

189.05 

0.175 

1.001 

0.0048 

0.00300 

0.21770 

490 

193.60 

0.121 

0.963 

0.0031 

0.00284 

0.23674 

510 

188.41 

0.103 

1.006 

0.0026 

0.00346 

0.34070 

555 

185.90 

0.076 

1.144 

0.0019 

0.00330 

0.35666 

670 

152.82 

0.137 

1.463 

0.0008 

0.00238 

0.29579 

765 

123.27 

0.040 

1.732 

0.0005 

0.00221 

0.32081 

865 

100.34 

0.010 

1.732 

0.0003 

0.00206 

0.34400 


The values for AT W (^) were taken from Baker and Smith 
(1982), and again spectrally weighted FWHM assuming a 
Gaussian response (Table 3). F C (A) was computed from 
Baker and Smith (1982): 


A C (A) = k c (X)[chl. a] exp 
-r 0.001 [chL a 


-KW logic 
2 


CVef J 


where [chi. a] is chlorophyll concentration, C re f is a refer- 
ence chlorophyll value (0.5), and k c ( A) and A:'( A) are spec- 
tral fit coefficients and were spectrally weighted over the 
SeaWiFS bands (Table 3). For this simulation, K g { A) was 
assumed to be zero. 

The backscattering coefficient b b (X) may be expressed 
using components similar to the A^(A) formulation (4): 


b b { A) = 0.5 bw(X) + b bc { A) (7) 

where bw{ A) is the total scattering coefficient for pure sea- 
water and b bc ( A) is the backscattering coefficient of phyto- 
plankton. The coefficient 0.5 is the backscattering-to-total 


scattering ratio for pure seawater. The backscattering co- 
efficient for plankton, fy> c (A), was determined from an em- 
pirical relationship developed by Gordon et al. (1988a) 

M A) = A(X)[chl. a] S(A) (8) 

where, the coefficients A(A) and B( A) were computed for 
SeaWiFS bands as shown in Table 3. The only remain- 
ing unknown for the model, then, is [chi. a], which was 
obtained from CZCS data. 

3.2 Orbit Model 

A simulated SeaWiFS orbit is required to obtain real- 
istic viewing and solar geometries, which are then used to 
produce the atmospheric contribution to the total radiance 
detected by the sensor. The method used here is to propa- 
gate a SeaWiFS orbit over a field of CZCS pigment values, 
use the values to compute L\ vat (A) for SeaWiFS bands, and 
use the simulated scan and tilt positions to determine the 
location of the pigment value in Earth coordinates (geo- 
centric latitude and longitude). Given this information, 
the spacecraft and solar zenith and azimuth angles can be 
calculated and then used to simulate the radiative prop- 
erties of an atmosphere. This section describes the orbit 
model employed, the equations to compute latitude and 
longitude, and the viewing and solar geometries. 

The orbital model used assumes a circular orbit. Thus, 
the only parameters required to propagate a simulated or- 
bit were inclination, altitude, and epoch (see Table 1 for 
SeaWiFS inclination and altitude). The epoch was chosen 
such that the first orbit crossed the equator on the de- 
scending node at the Greenwich Meridian at exactly noon 
local time. 

In such a simple orbit model, the location of succes- 
sive sub-satellite points (the point on the Earth where the 
projection of the satellite lies) is a simple function of the 
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previous sub-satellite point and the travel time between 
the two points. From spherical trigonometry, 

sin [#$(£ — t 0 )] = costf sin [# s (t 0 )] 

+ sin <5 cos [# s (£o)] cos(i') ^ 

where # s is the sub-satellite latitude at times t 0 and t-t 0 , 
6 is the great circle distance from # s (£ 0 ) to # s (£ — to), and 
i f is inclination i — 90° . The great circle distance is known 
from 

s = 360^ (10) 

where P is the nodal period. The longitude difference, 
Au 5 , is calculated similarly: 

. (ii) 

The new longitude, equals the old longitude u>o mi- 
nus the longitude difference Au> plus a correction for Earth 
rotation 

u = ujq — Aid + 0.004167 1 (12) 

w'here the coefficient 0.004167 is the rate of rotation (° s _1 ). 

3.3 Navigation of Pixels 

Given sub-satellite longitude and latitude, and the scan 
and tilt angles, the Earth position of any given pixel is cal- 
culable. First, one must calculate the great circle distance 
from the sub-satellite point to the pixel: 

£ = 6-e' s (13) 

where 0 is the spacecraft zenith angle, and 0' is the scan 
angle 8 S modified by tilt. Given 0 S , 0' is computed by 
a simple pitch rotation in the direction of the tilt. The 
spacecraft zenith angle, 0, is computed directly from 


A Wc — 


cos 


-1 


cos 6 — sin [# s (t — tp)] sin [# 5 (t 0 )] 


COS [# s (t - t 0 )] COS [# s (t 0 )] 


The pixel longitude is determined from Au; by either 
adding to, or subtracting from, the sub-satellite longitude, 
depending on which side of the scan the pixel is on; if it is 
west of the sub-satellite point a subtraction is performed. 

3.4 Viewing and Solar Geometries 

Once the pixel has been navigated to Earth coordi- 
nates, computation of viewing geometry (spacecraft zenith 
and azimuth angles) is straightforward. In fact, the space- 
craft zenith angle 0 has already been found (14). This 
leaves the computation of the spacecraft azimuth angle, 
which is defined as the angle from the pixel to the sub- 
satellite point, measured counterclockwise from true north 
(at the pixel) 


_ sin # s — cos 6 sin # 

# = — t (17 

sin 6 cos # 

where # is the spacecraft azimuth angle (the other terms 
are defined above). 

The determination of the solar geometry requires addi- 
tional information to the pixel location. Specifically, this 
additional information includes the solar declination lat- 
itude #<*, equation of time T e [both obtained using stan- 
dard methods, e.g., Iqbal (1983)], and the Greenwich Mean 
Time (GMT). GMT is determined by 

#gmt = *ec - + t € (18) 

where Hq mt is GMT in hours, £ec is the equator crossing 
time (in local time, i.e., noon) u e is the equator crossing 
longitude, and t € is the time difference to or from the most 
recent equator crossing to the present position (in hours). 

Given this information, plus the pixel location, the so- 
lar zenith angle 0 q is computed from 


0 


.-1 


R e + H s 
Re 


sin0; 


(14) 


cos 0o = sin # sin # d + cos # cos #<* cos 
where ft is the solar hour angle, computed from 


(19) 


where R € is the mean Earth radius (6371.2 km) and H s is Q = (H - 12)15 + w + T e . 

the spacecraft altitude (705km). 

The pixel latitude and longitude are determined from The solar azimuth angle $ 0 is computed by 
spherical trigonometry. For the latitude: 

# = sin” 1 [sin# 5 cos 6 + cos# 5 sin S cos 77] (15) = cos 


sin #<f - sin # cos 0 O 
cos # sin 0 O 


(20) 


( 21 ) 


where 77 is the bearing from the sub-satellite point to the 
pixel along the direction of motion of the satellite (either 
270° 4* i' for the left-hand side of the scan facing north or 
90° — i' for the right-hand side facing south). 

The longitude difference, Au;, from the sub-satellite 
point to the pixel is determined by 


The viewing geometries 0 and $ and solar geometries 
0o and $0 sire now known and attention may be focused 
on the addition of an atmosphere. 

3.5 Atmospheric Contributions 


Aid = cos 


cos 6 - sin # s sin# 
cos # s cos# 


This section focuses on the computation of realistic at- 
(20) mospheric contributions to the total radiance, based on 
representative atmospheric conditions and the simulated 
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SeaWiFS viewing and solar geometries. This means de- 
termining the Rayleigh scattering, ozone absorption, wa- 
ter vapor absorption, oxygen absorption, aerosol scatter- 
ing and absorption, and sun glint contributions. As a first 
step, absorption coefficients for the atmospheric gases must 
be determined and spectrally weighted over the SeaWiFS 
bands, as was done for the mean extraterrestrial irradiance. 
Absorption coefficients for ozone, -water vapor, and oxygen 
(the spectral absorption characteristics of which all three 
fall to some extent within some SeaWiFS bandwidths) are 
shown in Table 4, along with the Rayleigh optical thickness 
for the SeaWiFS bands. 

Table 4. Spectrally -weighted (FWHM) Rayleigh 
optical thickness (r r ), ozone absorption coefficient 
(a oz ), water vapor absorption coefficient (a wv ), and 
oxygen absorption coefficient (a ox ). Wavelengths, 

A, are in nm. The units for the absorption coeffi- 
cients are in cm -1 ; r r is dimensionless. 


A 

Tr 

&oz 

&WV 

^ox 

410 

0.3191 

0.0000 

0.0000 

0.0000 

443 

0.2364 

0.0027 

0.0000 

0.0000 

490 

0.1562 

0.0205 

0.0000 

0.0000 

510 

0.1326 

0.0382 

0.0000 

0.0000 

555 

0.0938 

0.0898 

0.0000 

0.0000 

670 

0.0437 

0,0463 

0.0006 

0.0000 

765 

0.0255 

0.0083 

0.0000 

6.9900 

865 

0.0155 

0.0000 

0.0008 

0.0000 


Following Gordon et al. (1983), the effects of the ab- 
sorbing gases in the atmosphere were determined by as- 
suming two trips of the extraterrestrial irradiance through 
the atmosphere. This was accomplished by 


The next step is to compute the Rayleigh scattering 
contribution L r to the total radiance received at the satel- 
lite. Since multiple scattering Rayleigh tables are not yet 
available for SeaWiFS bands, the tables for the CZCS 
bands are used (Gordon et al. 1988b). There are two co- 
incident bands for SeaWiFS and the CZCS (see Table 1): 
one at 443 nm and the other at 670 nm. The CZCS multi- 
ple scattering Rayleigh tables are used for these two bands. 
Values for non-coincident SeaWiFS bands are derived as- 
suming a A -4 principle. For bands 1-5 (412-555 nm), the 
procedure is as follows. Let A = 443 nm; then, 

x = /(A) A 4 . (25) 

Now, for A = 412 and 555 

/(A) = xA' 4 (26) 

from which the Rayleigh scattering is derived: 

Lr( A) = /(A)F'(A) (27) 

where /(A) is the Rayleigh intensity taken from the CZCS 
multiple scattering tables (Gordon et al. 1988b). 

The procedure is similar for bands 6-8 (670-865 nm), 
except 670 is substituted for A in (25), and 670 and 865 
replace 412 and 555 in (26). When this method is applied 
to the CZCS bands at 520 and 550 nm, a maximum error 
of 3.2% is obtained for 0 0 = 60 and 6 = 50. The error is 
less at all smaller angles tested. 

Next, the aerosol radiances, L a , must be determined. 
According to Gordon and Castano (1989) the aerosol ra- 
diance is determined by 


Ft = Fo T O2 (0 o) Toz(0) T wv (6q) T wv (0) T ox (0o) T ox (0) (22) 


where wavelength dependence has been dropped and T 
represents the transmittance of the gaseous component in 
subscript. The 0 O and 0 dependencies represent the two 
paths through the atmosphere: 6q is the path from the 
sun to the Earth and 0 is the path from the Earth to the 
satellite. The path lengths for each path are computed for 
water vapor and oxygen by (Kasten 1966): 


M 


1 

cos T) + 0.15(93.885 - r?)" 1 - 253 


(23) 


= i£g Tg Pa -Fp M (2g) 

47T 

where u> Q is the single scattering albedo of the aerosol, r a 
is the optical thickness, and p a is a factor to account for 
the probability of scattering to the spacecraft for three 
different paths from the sun (Gordon and Castano 1989). 
In (28) M has been substituted for Gordon and Castano’s 
(1989) l/cos0 to allow calculations on a curved Earth. 

The assumption is made that the aerosol obeys the 
Angstrom formulation 


where r\ is 0o or 0. 

The path length for ozone transmittance is slightly dif- 
ferent due to the high altitude of ozone concentrations 
(Paltridge and Platt 1976): 


■AF>2 — 


1.0035 

cos 2 7] -f 0.007 -0 5 


(24) 


These formulations avoid the assumption of a flat Earth. 


r a = 0 \~* (29) 

and the Gregg and Carder (1990) model may be used to 
obtain 3 and a, as a function of wind speed. Wind speeds 
were obtained from six years of data from the Fleet Numer- 
ical Oceanography Center (FNOC), made available from 
the NASA Climate Data System (NCDS). The mean wind 
speeds over these six years at 2.5 x 2.5° spatial resolution 
were used. 
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The aerosol scattering albedo, w a , was also computed Some of the scattered Rayleigh and aerosol components 
using Gregg and Carder (1990). The factor p a was com- continues forward and contributes to the total irradiance 
puted from Gordon and Castano (1989) at the surface. These are expressed separately as: 


Pa = 


4 7T 


[W) + [p(0o) + p(0)]P(0 + )] (30) 


t S — ^OZ ^WV t OX £ aa 


1 - 1 


0.95 


+ t}.- 5 F a (l -U 


(36) 


where p(0o) and p{6) are the Fresnel reflectances for solar 
and viewing geometries, respectively, and P is the aerosol 
scattering phase function for the forward scattering angles 
P(0 + ), i.e., scattering toward the surface and reflected off 
the ocean or reflected first and then scattered toward the 
spacecraft, and backward angles P(#“), i.e., backscattered 
from the atmosphere. 

P(#*) is determined using a Henyey- Greens tein func- 
tion for marine aerosols 


p(0±) = 


a(l ~9i) 

(1 + g\ — 2g\ cos®*) 1 - 5 
(l -q)(l -gp 
(l+^-2fl2Costf*)»-s 


(31) 


where a = 0.983, g\ — 0.82, and g*i =• —0.55 (Gordon and 
Castano 1989). As defined in Gordon and Castano (1989): 

cos#* — ±cos#ocos# — sin#o sin# cos(3> — <£o). (32) 

The next contribution to be accounted for is that of 
sun glint radiance diffusely transmitted to the spacecraft 
(TL g ). Sun glint at the surface for a given viewing and 
solar geometry may be expressed by (Cox and Munk 1954): 

L = Pn To(A, # 0 ) Pq(A) p w (0, g, #q, gp, W) 

9 4cos#cos 4 #^ 

where p t v (0, <£,#o, $o) is the probability of seeing sun glit- 
ter in the direction #, given the sun in position # 0,^0 as 
a function of wind speed (W\ and 9 N is the angle with re- 
spect to nadir of the sea surface slopes to produce a reflec- 
tion angle to the spacecraft (Viollier et al. 1980). The term 
To(A,#o) represents the total downward transmittance of 
irradiance as a function of #0 and the absorbing gases in 
the atmosphere. It may be expressed as the sum of the 
direct and diffuse transmittances (again dropping wave- 
length dependence) 


w here £ aa is the aerosol transmittance after absorption 

*aa = e~ {1 ~^ a)TaM (37) 

and is analogously the transmittance after scattering 
(Justus and Paris 1985) 

(38) 

In (36)-(38), F a is the forward scattering probability of 
the aerosol, taken from Bird and Riordan (1986), w a is the 
single scattering albedo of the aerosol, and r a is the aerosol 
optical thickness. 

If L g ( A) is the sun glint radiance at the surface, TL g ( A) 
is that received by the sensor, where T is the total trans- 
mittance from the Earth to the satellite. It is defined in 
(34)-(38) with # substituted for # 0 . 

Finally, the normalized water-leaving radiance Lwn 
must be converted into the water-leaving radiance mea- 
sured at the satellite TL\y. This is expressed by Gordon 
(1990) as 


TL\y = L\yn ( I ~ ^)To cos#o. (39) 

The weighted direct plus diffuse reflectance, p, is computed 
as: 

P = PnW d + p N W s (40) 

where W d equals the direct irradiance divided by the total 
irradiance at the surface, and W s is the diffuse irradiance 
divided by the total. These components are computed as a 
function of the atmospheric constituents and solar zenith 
angle using the model of Gregg and Carder (1990), as are 
the direct and diffuse reflectances, which are functions of 
the wind speed. 

The simulated total radiance at the spacecraft, given 
in (1) is calculated by summing (27), (28), (33), and (39). 

3.6 Atmospheric Conditions 


r 0 (#o) = td(0o) 4- t s (9 o) (34) 

where the subscripts d and s represent direct and diffuse 
components, respectively. The transmittance t d is simply 
the transmittance after absorption by the gaseous com- 
ponents of the atmosphere, scattering and absorption by 
aerosols, and scattering by Rayleigh (dropping the angular 
dependence) : 

td t r t aa^oz^wv^ox* (35) 


There are five attenuating sources in the atmosphere 
that lie within SeaWiFS bands: Rayleigh scattering, ozone 
absorption, water vapor absorption, oxygen absorption, 
and aerosol scattering and absorption. Each w^as taken 
for the simulated data set at mean values. Rayleigh scat- 
tering and oxygen absorption are functions of atmospheric 
pressure for which a standard value of 1,013.25 mb was 
used. Ozone absorption was taken to be 340 Dobson units 
(DU) and absorption due to water vapor was assumed to 
be 1.5 cm. 
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Aerosols may be calculated as a function of visibility, 
a commonly reported meteorological parameter. Visibility 
was set to 15 km for the simulation, which corresponds 
to a mildly turbid atmosphere. Relative humidity (which 
determines, in part, the size distribution of the aerosols) 
was set at 80%. The 24-hour mean wind speed required for 
the aerosol model used by Gregg and Carder (1990), was 
estimated at 4.75 m s" 1 , the global mean from 6 years of 
FNOC data, and the current wind speed, as noted before, 
was determined as the mean over 2.5° x2.5° grid cells from 
the 6- year FNOC data set. 


Table 5. Lner and saturation radiances (L sat ) for 
SeaWiFS (units are mW cm -2 /im' 1 sr -1 ). 


Band 

A [nm] 

Lner 

Tsat 

1 

412 

0.0182 

13.63 

2 

443 

0.0125 

13.25 

3 

490 

0.0098 

10.50 

4 

510 

0.0088 

9.08 

5 

555 

0.0077 

7.44 

6 

670 

0.0056 

4.20 

7 

765 

0.0035 

3.00 

8 

865 

0.0023 

2.13 


3.7 Ten-bit Words and Data Structures 4 APPLICATION OF CZCS DATA 


The final step in developing a realistic SeaWiFS sim- 
ulated data set is to convert the science data into 10- bit 
words and to structure the data fields according to the 
expected format for the downlink data. Not all of the 
telemetry fields have been decided upon, but the purpose 
here is to create as reasonable a data structure as possible 
to allow data system testing. Further developments can be 
incorporated as they are made. 

Fig. 1 show f s the preliminary stored data structures for 
LAC and GAC data. Depicted are the structures for a 
major frame of GAC and a major frame of LAC, each 
consisting of 3 minor frames. A major frame of GAC thus 
consists of 15 GAC scan lines, w r hile a major frame of LAC 
consists of 3 LAC scan lines. The number of bytes in each 
minor frame is also denoted. 

A single LAC scan line of SeaWiFS data contains 1,288 
pixels, comprising 1,285 science pixels, a start synchroniza- 
tion pixel (also a 10-bit w'ord) preceding the first pixel, a 
stop synchronization pixel at the end of the scan line, fol- 
lowed by a time delay integration (TDI) pixel. A GAC 
scan likewise contains 251 pixels, w*ith 248 actual data pix- 
els and three others in the same format as the LAC data. 
This representation is included in the simulated data set. 

SeaWiFS simulated total radiance data were packed 
into 10-bit words according to 

DC, 0(A) = (41) 

$(A) 


where DCio is the digital counts at 10-bit digitization, 
Lner(^) is the noise equivalent radiance (NER), or the 
minimum detectable radiance, and s is the slope for the 
range 0-1,023, wdiich is applicable for 10-bit digitization, 
given by 


s(A) 


Tsat (A) - LnerM 
1,023 


(42) 


where L sat ( A) is the saturation radiance for the sensor as a 
function of wavelength. Tner(^) and T sat (A) are showm in 
Table 5. This provides a realistic data set corresponding 
to the format stored on board the data recorder, and thus 
may be considered a level-0 data set. 


The SeaWiFS GAC simulated data set w'as derived 
from the global CZCS data set (Feldman et al. 1989). This 
level-3 (remapped and gridded) data set w*as in dimensions 
of 2,048 x 1,024 (0.176° longitude by 0.176° latitude), and 
w'as cloud- and land-masked. The simulated SeaW r iFS or- 
bits passed through this CZCS pigment field in descending 
node (from north to south) for approximately 6.8 orbits to 
simulate a single data recorder dump of data. Start time 
for the simulation w'as the first contact at Wallops Flight 
Facility (WFF), the dowmlink station for SeaWiFS GAC 
data. Succeeding start times, as w'ell as all stop times, w'ere 
derived in previous simulations to maximize GAC cover- 
age. The ultimate limit of the sensor was assumed to be 2 
NER values, and the recorder w'as turned on whenever this 
limit w'as reached, subject to land features and ice cover. 
An 80% ice concentration w T as assumed as the limit for re- 
mote sensing purposes. Ice concentrations used w'ere the 
minimum of six years of International Satellite Cloud Cli- 
matology Project (ISCCP) data, made available by NCDS. 
The start and stop times, in minutes travelled along the 
orbital track from equator crossing at ascending node, are 
given in Table 6. 

The tilt was left at nadir pointing (0°) for this simu- 
lation to test the response of the data system to excessive 
sun glint. Normally, the sensor wall be tilted aft w'hen 
approaching the solar declination latitude from the north, 
and turned forward for the remainder of the descending 
node track; however, this strategy' has yet to be finalized. 

For the GAC simulation, no attempt was made to in- 
clude clouds, and land features w'ere masked (set to zero) 
as they occurred in the CZCS global data set. 

The LAC simulation was initiated with the same global 
CZCS pigment data, using the same orbit trajectories, as 
will be the case for the real spacecraft. The sensor LAC 
recorder was turned on for three separate occasions, com- 
prising one 3-minute scene and two 2-minute scenes, for 
a total of 7 minutes of LAC record time. This is about 
the correct amount of time left on the recorder after GAC 
data recording and some calibration activities. The start 
and stop times for the LAC recorder, and the associated 
simulated data set, are provided in Table 7. The difference 
between LAC and GAC data is the amount of data in a 
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Fig. 1. Structure of the GAG and LAC data streams as stored on board. S/C indicates spacecraft, TLM means telemetry, and SOU means state of health. 
The number of words in each field is represented above the field, and the number of bytes in each field is represented underneath. The entries in slanted 
typeface are in an 8-bit telemetry format and the plain typeface entries are in a 10-bit format. (Courtesy of Steven Jackson, OSC.) 
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Table 6. Start and stop times (measured in minutes travelled along the orbit track from the equator crossing 
on ascending node), total time, start and stop latitudes, and solar zenith angles for the orbital simulations for 
the SeaWiFS simulated data set. The first contact with WFF is made on orbit 4 at 33.8743 min. and the second 
contact is on orbit 11 at night. Thus, the simulated data begins at 33.8743 min. and ends at the stop time of 
orbit number 10. 


Day 

Orbit 

Number 

Times 

Latitude 

Solar Zenith 

Start 

Stop 

Total 

Start 

Stop 

Start 

Stop 

80 

1 

27.0571 

71.3008 

44.243 

78.285 

-76.530 

81.396 

79.507 

80 

2 

26.9876 

69.2139 

42.226 

78.460 

-69.996 

81.648 

71.923 

80 

3 

26.9043 

69.9826 

43.078 

78.670 

-72.499 

81.954 

74.716 

80 

4 

26.9154 

70.7818 

43.866 

78.642 

-74.995 

81.910 

77.621 

80 

5 

27.7986 

70.5321 

42.733 

76.235 

-74.230 

78.701 

76.713 

80 

6 

26.8932 

71.6449 

44.751 

78.698 

-77.496 

81.991 

80.757 

80 

7 

26.9099 

71.6449 

44.735 

78.656 

-77.496 

81.931 

80.757 

80 

8 

26.8932 

69.2139 

42.320 

78.698 

-69.996 

81.991 

71.923 

80 

9 

30.3202 

68.4647 

38.144 

68.207 

-67.492 

69.537 

69.200 

80 

10 

28.7511 

68.4647 

39.713 

73.336 

-67.492 

75.239 

69.200 

80 

11 

27.8374 

68.4647 

40.627 

76.121 

-67.492 

78.560 

69.200 

80 

12 

26.8932 

68.4647 

41.571 

78.698 

-67.492 

81.991 

69.200 

80 

13 

26.9876 

| 69.2139 

42.226 

78.460 

-69.996 

81.648 

71.923 

80 

14 

27.0571 

69.2139 

42.156 

78.282 

-69.996 

81.396 

71.923 


LAC scan (Table 1), which is 1,285 pixels, with no gaps 
along-scan. 

Table 7. Start and stop times (in minutes trav- 
elled along the orbit track from the equator crossing 
on ascending node) orbit numbers, and correspond- 
ing latitudes for the LAC simulation. 


Orbit 

Number 

Times 

Latitude 

Start 

Stop 

Total 

Start 

Stop 

4 

54.6 

57.6 

3.0 

-18.35 

-29.13 

6 

64.3 

66.3 

2.0 

-53.03 

-60.06 

6 

68.6 

70.6 

2.0 

-67.95 

-74.44 


5. RESULTS AND DISCUSSION 

Imagery of the simulated SeaWiFS total radiances L t (A) 
for all 8 bands is shown in Figure 2 for GAC data. The 
start of the simulation was at the first visibility at WFF, 
which can be noted by the truncated orbit at the east- 
ern edge of the image. Recall that for this simulation the 
sensor remained untilted, resulting in substantial sun glint 
influence near the equator (which is the solar declination 
for the vernal equinox, the time of the simulation). Tilting 
the sensor aft of the velocity vector when approaching the 
solar declination, and then forward when leaving would 
decrease the sun glint contribution to the total radiance 
substantially. 

Note the decreases in L t ( A) toward the poles, due to 
a reduction in the flux of solar irradiance into the oceans, 
and the increases in L t { A) at the scan edges, due to in- 
creased atmospheric contribution. The bands show pro- 
gressively less total radiance from 412-865 nm. However, 
band 7 (765 nm) is the lowest because of large oxygen ab- 


sorption (see Table 4). This absorption reduces the amount 
of saturation in the sun glint regions. Little evidence of 
water-leaving radiances is evident in most of the images. 
This is because the atmosphere dominates the total ra- 
diance over the oceans. Some evidence may be seen in 
the image of band 2 (443 nm), where the contribution of 
water-leaving radiance to the total is greatest under low 
chlorophyll concentrations. 

The LAC recording areas are also depicted in Figs. 2-4. 
Note that they correspond with portions of the GAC or- 
bits, but have a wider swath. Note also that the extra data 
across-track have very large L t { A) values, due to increased 
atmospheric contribution as a function of the greater LAC 
spacecraft zenith angle at the edges. 

This first version of the SeaWiFS simulated data set 
contains realistic total radiances, in a representative data 
record scenario: approximately 6.8 orbits of GAC data be- 
ginning at the first visibility of the dowmlink station at 
WFF, and 8 minutes of LAC record time, interspersed in 
the SeaW r iFS orbits. Future improvements should incorpo- 
rate representative values in the telemetry fields, account- 
ing for clouds and land radiances. Also, a tilt strategy 
to minimize sun glint should be included. However, this 
version facilitates early prototyping activities by the Sea- 
WiFS data system, providing a realistic data set to assess 
performance. 

6- DATA AVAILABILITY 

Version 1 of the simulated SeaWiFS data set is avail- 
able from the authors on 9-track tape at 6,250 bits per inch 
(bpi). The tape is written using the VAX/VMS BACKUP 
command, on the assumption that most users have access 
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Fig. 2. Simulated total radiances for GAC data for SeaWiFS bands 1-3. 
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Band 4 



Band 5 



Band 6 

Fig. 3. Simulated total radiances for GAC data for SeaWiFS bands 4-6. 
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Band 2 LAC 

Fig. 4. Simulated total radiances for GAC data for SeaWiFS bands 7-8, and LAC data for band 2. 
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to Digital Equipment Corporation (DEC) hardware. FOR- 
TRAN programs to read the data and strip out telemetry 
fields and science data are provided in the Appendices. 

Appendices 


A. Description 

B. Program RGAC.FOR 

C. Program RLAC . FOR 

D. Program EXPAND. FOR 

E. Program B2W. FOR 

Appendix A 

Description 

The SeaWiFS simulated data sets and programs to read 
them are provided on a 6,250 bpi 9-track tape. The tape was 
generated using the VMS (version 5.4-3) BACKUP command and 
the following command sequence: 

MOUNT/FOREIGN device 

BACKUP/LOG device : simdat * . */label=simdat/block=32256 

DISMOUNT device 

where device is the tape drive device name. (Note: 105 MBytes 
of free disk space are required to read the tape in its entirety.) 

When the tape is loaded onto disk, 6 files are put in the 
current directory. GAC . DAT and LAC . DAT are the simulated data 
sets for GAC and LAC data, respectively. Listings of the 
source codes to read these data sets are provided in the follow- 
ing. This includes the main programs RGAC.FOR and RLAC. FOR 
and the two supporting subroutines, EXPAND. FOR and B2W.F0R. 
EXPAND. FOR is used to extract the minor frame data segments. 
B2W.F0R is used to unpack a 5- byte field into four 2- byte inte- 
gers. 


Appendix B 
Program RGAC.FOR 


program rgac 
c 

c Program to read and sort simulated SeaWiFS GAC data, 
c by: Frank Chen, GSC. 
c 

implicit none 
character c( 13864) 
integer*2 si (12) 
character s2(775) 
integer*2 s3( 10458) 
character s4(2) 
c 

c Minor frame data structures 
c 

integer*2 FrameSync(6) ,SCID(2) ,TimeTag(4) 
character SCS0HC775) 

integer* 2 InstTLM (44) , Callable (44) ,Spare(188) 
integer*2 GAC(8,252,5) 
integer*2 PAD 1(2) , AuxSync (100) 
character PAD2(2) 
c 

c Science data array 
c 

integer*2 ilt(8,248) 
c 

c Total number of scan lines 
c 

integer*4 nscan 
c 

c Minor frame number (1 ,2,3) 


c 

c integer*4 MinorFrameNo,msec_of _day , spacecraf t_id 
integer*2 day_of_year 
integer*4 i,j,k 
c 

c Orbit number, xlons, ylats, and gmt time (hr) for 5 GAC 
c scan lines within one minor frame. These fields are 
c temporary, 
c 

integer*4 iorb(5) 

real*4 xlons(5) ,ylats(5) ,gmt (5) 
c 

c Map the above fields within SCS0H 
c 

equivalence (iorb(l) , SCS0H( 197) ) , (xlons ( 1) ,SCS0H (217) ) 
equivalence (ylats( 1) ,SCS0H(237) ) , (gmt (1) ,SCS0H(257) ) 
c 

c Map fields within minor frame to the three segments 
c 

equivalence (FrameSync(l) , si (1) ) , (SCID(l) ,sl(7)) 
equivalence (TimeTag(l) ,sl (9)) , (SCS0H( 1) , s2(l) ) 
equivalence (InstTLM(l) ,s3(l)) , (CalTable(l) ,s3(45) ) 
equivalence (Spared) ,s3(89)) , (GAC(1 , 1 , 1) ,s3(277)) 
equivalence (PAD1(1) ,s3( 10357)) , (AuxSync (1) ,s3( 10359)) 
equivalence (PAD2(1) ,s4(l)) 
c 

100 f ormat (i5 ,3f 10 .4) 

200 format (8i5) 
nscan - 0 
c 

open (4 ,f ile= ’gac .dat ’ , status=’old J ,fonn=’ unformatted’) 

20 continue 


Read in one minor frame 


read (4 , end=999) c 
c 

c Convert all 10-bit-word to 16-bit-word 
c 

call expand(c,sl,s2,s3,s4) 
c 

c Calculate s/c id, minor fram number, day of year, 
c and msec of day 
c 

spacecraf t_id = M0D(SCID(l)/8 , 16) 

MinorFrameNo - M0D(SCID(1)/128,4) 
day_of _year = TimeTag(l)/2 
msec_of_day = TimeTag(2)*1024 .0*1024 . 0 + 

* TimeTag(3)*1024. 0 + TimeTag(4) 


do 

c 

c 

c 

c 

c 

c 

c 

c 


i=l , 5 

GAC(k, 

GAC(k, 

GAC(k, 


1 . 1) ,k=l ,8 

2.1) ,k=l,8 

3.1) ,k=l,8 
GAC (k, 252, i) ,k-l ,8 


80 bits GAC GainTDI 
80 bits GAC StartSync 
80 bits GAC DarkRestore 
80 bits GAC StopSync 


Check if the GAC data are valid 


if ( (GAC(1 ,3, i) .eq.O) .and. 
(GAC(3,3,i) .eq.O) .and. 

(GACC5 ,3, i) .eq.O) .and. 

(GAC(7 ,3, i) .eq.O) .and. 
write(6,*)’ no more valid data’ 
goto 999 
end if 

nscan = nscan + 1 
do j=4, 251 


<GAC(2 ,3 , i) . eq.O) 
(GACC4 ,3 ,i) . eq.O) 
(GAC(6,3,i) .eq.O) 


. and. 
. and. 
. and. 


(GACC8 ,3 , i) . eq.O) ) then 


do k=l , 8 

ilt (k , j-3) s =GAC(k, j , i) 
end do 
end do 
end do 
goto 20 
999 continue 
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close(4) 

write(6,*) 'Total no. of scan lines nscan 
c 

end 

Appendix C 
Program RLAC . FOR 

program rlac 
c 

c Program to read and sort SeaWiFS simulated LAC data, 
c by: Frank Chen, GSC. 
c 

implicit none 
character c (13864) 
integer*2 sl(12) 
character s2(775) 
integer*2 s3( 10458) 
character s4(2) 
c 

c Minor frame data structures 
c 

integer*2 FrameSync(6) ,SCID(2) f TimeTag(4) 
character SCS0HC775) 
integer*2 InstTLM (44) , LAC (8 , 1289) 
integer*2 PAD1 (2) , AuxSync (100) 
character PAD2(2) 
c 

c History SCSOH and 0X25 in minorframe 2,3 
c 

character HistS0H(775) , 0X25(775) 
c 

c Ancilliary TLM and Calibration Table in minor 
c frame 2,3 
c 

integer*2 AncTLM(44) ,CalTable(44) 
c 

c Science data array 
c 

integer*2 ilt (8 , 1285) 
c 

c Total number of scan lines 
c 

integer*4 nscan 
c 

c Minor frame number (1 ,2,3) 
c 

integer*4 MinorFrameNo, msec_of _day , spacecraf t_id 
integer*2 day_of_year 
integer*4 i,j,k 
c 

c Orbit number ,xlons,yl at s , and gmt time (hr) for LAC 
c scan line within one minor frame. These fields 
c are temporary . 
c 

integer*4 iorb 

real*4 xlons,ylats ,gmt 

c 

c Map the above fields within SCSOH 
c 

equivalence ( iorb, SCSOH ( 197) ) , (xlons,SCS0H(217) ) 
equivalence (ylats, SCSOH (237)), (gmt , SCSOH (257) ) 
c 

c Map fields within minor frame to the three segments 
c 

equivalence (FrameSync(l) ,sl(l)) , (SCID(l) ,sl (7) ) 
equivalence (TimeTag(l) ,sl(9)) , (SCS0H(1) ,s2(l) ) 
equivalence (InstTLM(l) ,s3(l)) , (LAC(1,1) ,s3(45)) 
equivalence (PAD1(1) ,s3(10357)) , (AuxSync ( 1) , s3( 10359) ) 
equivalence (PAD2(1) t s4(l)> , (HistSOH(l) ,SCS0H(1)) 
equivalence (0X25(1) ,SCS0H(1)) , (AncTLM(l) ,InstTLM(l)) 
equivalence (Callable (1) ,InstTLM(l)) 
c 


100 format (i5,3f 10.4) 

200 format(8i5) 
nscan = 0 
c 

open (4,file= 5 lac .dat ' »status= 5 old 5 , form= 'unformatted 5 ) 
20 continue 
c 

c Read in one minor frame 
c 

read(4,end=999)c 

c 

c Convert all 10-bit-word to 16-bit-word 
c 

call expand(c,sl,s2,s3,s4) 
c 

c Calculate s/c id, minor fram number, day of year, 
c and msec of day 
c 

spacecraft.id = M0D(SCID(1) /8, 16) 

MinorFrameNo = M0D(SCID(1)/128,4) 

day_of _year = TimeTag(l)/2 

msec_of _day = TimeTag (2) *1024.0*1024 .0 + 

* TimeTag(3) *1024.0 + TimeTag(4) 

c 

c write(6,*) 'day: 5 ,day_of _year , 5 scid: 5 , 
c * spacecraft_id, 

c * 5 f rame : 5 , MinorFrameNo 

c 

c LAC(k, 1) ,k-l ,8 80 bits LAC GainTDI 

c LAC(k, 2),k=l,8 80 bits LAC StartSync 

c LAC(k, 3) ,k=l ,8 80 bits LAC DarkRestore 

c LAC (k, 1289) ,k=l,8 <80 bits LAC StopSync 
c 

nscan = nscan + 1 
do j=4 , 1288 
do k=l,8 

i It (k , j -3) =LAC (k , j ) 
end do 
end do 
goto 20 
999 continue 
close(4) 

write (6,*) 'Total no. of scam, lines nscan 
end 

Appendix D 

Program EXPAND. FOR 

c This subroutine converts the input buffer (c) into 
c four minor frame segments: si, s2, s3, and s4. 
c SI consists of 12 10-bit words containing the frame 
c sync, the spacecraft id, and the time tag. S2 
c is a the S/C S0H array. S3 is the remainder 10458 
c 10-bit words without the last padding (1.5 bytes), 
c S4 is the last padding, 
c 

subroutine expand(c,sl,s2,s3,s4) 

implicit none 

character c( 13864) 

integer*2 si (12) 

character s2(775) 

integer*2 s3( 10458) 

character s4(2) 

character b(5) 

integer*2 w(4) 

integer*4 i,j,k 

c 

c Process first 1*2 segment 
c 

do i=0,ll,4 
k = i*5/4 
do j=l,5 
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b(j) = c (k+j ) 
end do 

call b2v(b,v) 
do j=l,4 

sl<i+j) = w(j) 

end do 
end do 
c 

c Process second C*1 segment 
c 

k = 12*5/4 
do i=l ,775 

s2(i) = c(k+i) 
end do 
c 

c Process third 1*2 segment 
c 

do i=0, 10455,4 

k - (i+12)*5/4 + 775 
do j-1 ,5 

b(j) = c(k+j) 
end do 

call b2v(b,w) 
do j=l,4 

s3(i+j) = w(j) 
end do 
end do 
k = 13860 
do j-1, 3 

b(j) = c(k+j) 
end do 

b(4) = char Ci char (c(k+4)) .and. 240) 
b(5) * char(0) 
call b2v(b,v) 
do j*l,3 

s3( 10456+ j ) = w(j) 
end do 

s4(l) = char (MOD (ichar(c( 13863) ) , 16) *16 + 

* ichar(c(13864))/16) 

s4(2) = char (MOD (i char (c (13864) ) , 16) *16) 

return 

end 

Appendix E 

Program B2W.F0R 


c Convert packed 5 bytes character (b) into four 
c two-byte- integer (v) each contains 10-bit information, 
c M0D(n,2**m): modulus , n 1 by . extract lowest 

c ’m’ bits from the number 'n*. 

c n * (2**m) : shift left the number ’el’ by 3 m > bits, 
c n / (2**m) : shift right the number > n > by 5 m* bits, 
c 

subroutine b2u(b,v) 
implicit none 
character b(5) 
integer*2 w(4) 
integer*4 c(5) 
integer*4 i 

integer*4 pO /!/, p2 /4/, p4 /16/, p6 /64/, 
integer*4 p8 /256/ 
c 

do i*l,5 

c (i) = ICHAR(b(i)) 
end do 

v(l) * M0D(c(l) ,p8)*p2 + c(2)/p6 

w(2) = M0D(c(2) ,p€)*p4 + c(3)/p4 

w(3) = M0D(c(3) ,p4)*p6 + c(4)/p2 

w(4) = M0D(c(4) ,p2)*p8 + c(5)/p0 

return 
end 


bpi bits per inch 

CZCS Coastal Zone Color Scanner 

DEC Digital Equipment Corporation 
DU Dobson Units 

FNOC Fleet Numerical Oceanography Center 
FORTRAN Formula Translation (computer language) 

FWHM Full-Width Half-Maximum 

GAC Global Area Coverage 
GMT Greenwich Mean Time 

IFOV Instantaneous Field-of-View 
ISCCP International Satellite Cloud Climatology Project 

LAC Local Area Coverage 
Level-0 Raw data. 

Level-3 G ridded and averaged derived products. 

MB Megabytes 
MF Major Frame 
mF Minor Frame 

NCDS NASA Climate Data System 
NER Noise Equivalent Radiance 

OSC Orbital Sciences Corporation 

SeaWiFS Sea- viewing Wide Field-of-view Sensor 
S/C Spacecraft 
SOH State of Health 

TDI Time Delay Integration 
TLM Telemetry 

VAX Virtual Address Extension 
VMS Virtual Memory System 

WFF Wallops Flight Facility 

Symbols 

a A constant equal to 0.983. 
a ox Coefficient for oxygen absorption. 
a oz Coefficient for ozone absorption. 
a wv Coefficient for water vapor absorption. 

A (A) Coefficient for calculating bb( A). 

bb (A) Spectral backscattering coefficient. 
bbc( A) Spectral backscattering coefficient for phytoplank- 
ton. 

b w (X) Total scattering coefficient for pure seawater. 

B( A) Coefficient for calculating bb( A). 

Gref Reference chlorophyll value (0.5). 

[chi. a] Chlorophyll concentration. 

D Sequential day of the year. 

DCio Digital counts at 10-bit digitization. 

Fo(A) Mean extraterrestrial irradiance. 

Fo Extraterrestrial irradiance corrected for Earth- 
sun distance. 

Fo Extraterrestrial irradiance corrected for the at- 
mosphere. 

F a Forward scattering probability of the aerosol. 
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g\ A constant equal to 0.82. 
go A constant equal to —0.55. 

Hq m t GMT in hours. 

H s Altitude of the spacecraft (for SeaStar 705 km). 

i Inclination angle. 
x Inclination angle minus 90°. 

I Rayleigh intensity. 

k c ( A) Spectral fit coefficient weighted over the Sea- 
WiFS bands. 

Ki (A) Spectral fit coefficient weighted over the Sea- 
WiFS bands. 

K( A) Spectral attenuation coefficient. 

K c { A) Attenuation coefficients for phytoplankton. 

Ke{ A) Attenuation coefficient downwelled irradiance. 
K g ( A) Attenuation coefficients for Gelbstoff. 

Ki(z,X) Attenuation coefficient upwelled radiance . 

K w ( A) Attenuation coefficients for pure seawater. 

L a { A) Aerosol radiance. 

L g { A) Sun glint radiance. 

TnerA) Noise equivalent radiance. 

L r (A) Rayleigh radiance. 

Tsat(A) Saturation radiance for the sensor. 

L t ( A) Total radiance at the sensor. 

L u (z, A) Upwelled spectral rad iance. 

Lw{ A) Water-leaving radiance. 

Lwn{ A) Normalized water-leaving radiances leaving an 
ocean of given optical properties assuming no 
atmosphere and the sun directly overhead. 

M Path length through the atmosphere. 

A/ oz Path length for ozone transmittance. 

rt Index of refraction. 

Q Irradiance- to- radiance ratio (equals n for totally 
diffuse radiance). 

p a A factor to account for the probability of scat- 
tering to the spacecraft for three different paths 
from the sun. 

p w The probability of seeing sun glitter in the di- 
rection given the sun in position 0 o, $o as a 
function of wind speed ( W ). 

P Nodal period. 

P a Probability of scattering to the spacecraft. 
P{0 + ) Phase function for forward scattering. 

P(0“) Phase function for backward scattering. 

r Water-air reflectance for totally diffuse irradi- 
ance. 

R{ 0” , A) Irradiance reflectance just below the sea surface. 
R e Mean Earth radius (6371.2 km). 

s(A) Slope for the range 0-1,023. 

t Time variable. 
to Initial time. 

£ aa Aerosol transmittance after absorption, 
tas Aerosol transmittance after scattering. 
td Direct component of transmittance after absorp- 
tion by the gaseous components of the atmo- 
sphere, scattering and absorption by aerosols, 
and scattering by Rayleigh. 


t e Time difference in hours between present posi- 
tion and most recent equator crossing. 
tEC Equator crossing time, 
toz Transmittance after absorption by ozone. 
t r Transmittance after Rayleigh scattering. 
t s Diffuse component of transmittance after absorp- 
tion by the gaseous components of the atmo- 
sphere, scattering and absorption by aerosols, 
and scattering by Rayleigh. 
tw Transmittance after absorption by water vapor. 

T(A,0) Total transmittance (direct plus diffuse) from 
the ocean through the atmosphere to the space- 
craft along the path determined by the space- 
craft zenith angle 0. 

T(\,6)Lw(X) The water-leaving radiance transmitted to the 
spacecraft. 

To(A,#o) Total downward transmittance of irradiance. 

T e Equation of time. 

T 0 x Transmittance of oxygen (O 2 ). 

Toz Transmittance of ozone (O 3 ). 

T S (A) Transmittance through the surface. 

T wv Transmittance of water vapor (H 2 O). 

W Wind speed. 

Wd Direct irradiance divided by the total irradiance 
at the surface. 

W s Diffuse irradiance divided by the total irradi- 
ance. 

a The power constant in the Angstrom formula- 
tion. 

(3 A constant in the Angstrom formulation. 

6 Great circle distance from ty s (to) to ^ 3 {t — to). 
A£ Time difference. 

Au; The longitude difference from the sub-satellite 
point to the pixel. 

Ac j s Longitude difference. 

7 } Bearing from the sub-satellite point to the pixel 
along the direction of motion of the satellite. 

A Wavelength of light. 

'I' Pixel latitude, 

'I'd Solar declination latitude. 

Sub-satellite latitude as a function of time. 

Spacecraft azimuth angle. 

4>o Solar azimuth angle. 

p Weighted direct plus diffuse reflectance. 
p{6) Fresnel reflectance for viewing geometry. 
p(0 o ) Fresnel reflectance for solar geometry. 

p n Sea surface reflectance for direct irradiance at 
normal incidence for a flat sea. 
p;v Reflectance for diffuse irradiance. 

0 Spacecraft zenith angle. 

60 Solar zenith angle. 

Bn The angle with respect to nadir that the sea sur- 
face slopes to produce a reflection single to the 
spacecraft. 

0 9 Scan angle of sensor. 

B' s Scan angle of sensor adjusted for tilt. 
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r a Aerosol optical thickness. 
r r Rayleigh optical thickness. 
r s ( A) Spectral solar atmospheric transmission. 

u; Longitude variable, 
too Old longitude value. 

Lt/’a Single scattering albedo of the aerosol. 
u; e Equator crossing longitude. 
lu s Longitude variable. 

Q Solar hour angle. 
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